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ABSTRACT 

Existing models invoking AGN activty to resolve the cooling flow conundrum in galaxy 
ly^ ' clusters focus exclusively on the role of the central galaxy. Such models require fine- 

tuning of highly uncertain microscopic transport properties to distribute the thermal 
thermal over the entire cluster cooling core. We propose that the ICM is instead heated 
' by multiple, spatially distributed AGNs. The central regions of galaxy clusters are rich 

I in spheroidal systems, all of which are thought to host black holes and could participate 

. in the heating of the ICM via AGN activity of varying strenghts. And they do. There 

is mounting observational evidence for multiple AGNs in cluster environments. Active 
AGNs drive bubbles into the ICM. We identify three distinct interactions between 
the bubble and the ICM: (1) Upon injection, the bubbles expand rapidly in situ to 
' reach pressure equilibrium with their surroundings, generating shocks and waves whose 

, dissipation is the principal source of ICM heating. (2) Once inflated, the bubbles rise 

buoyantly at rate determined by a balance with the viscous drag force, which itself 
results in some additional heating. (3) Rising bubbles expand and compress their 
surroundings. This process is adiabatic and does not contribute to any additional 
Q ■ heating; rather, the increased ICM density due to compression enhances cooling. Our 

|H I model sidesteps the "transport" issue by relying on the spatially distributed galaxies to 

C/3 i heat the cluster core. We include self-regulation in our model by linking AGN activity 

' in a galaxy to cooling characteristics of the surrounding ICM. We use a spherically 

symmetric one-dimensional hydrodynamical code to carry out a preliminary study 
illustrating the eflicacy of the model. Our self-regulating scenario predicts that there 
should be enhanced AGN activity of galaxies inside the cooling regions compared to 
galaxies in the outer parts of the cluster. This prediction remains to be confirmed or 
refuted by observations. 

Key words: cosmology: theory, observation , dark matter, large-scale structure of 
the Universe — gravitation 
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1 INTRODUCTION 

The observed properties of the x-ray emitting gas in clus- 
ters of galaxies present several challenging puzzles for which 
satisfactory explanations are yet to be found. The absence 
of cold gas (with temperatures below ~ 1 keV) is one of 
those puzzles (e.g. Peterson et. al. 2001). Because of radia- 
tive losses, one would expect the temperature profiles to fall 
below a keV or so in the central regions where cooling is most 
efficient. Yet in none of the observed clusters does the tem- 
perature drop to the level dictated by cooling alone. Hence 



efficient heating mechanisms must operate without fail at 
the cores of all cooling clusters. 

The most popular mechanism for suppressing cooling is 
energy released by an AGN harboured by the central cluster 
galaxy (cf. Quilis et. al. 2001, Babul et. al. 2002, Kaiser & 
Binney 2003, Dalla Vecchia et. al. 2004, Omma et. al. 200, 
Roychowdhury et. al. 2004, Ruszkowski et. al. 2004, Voit 
& Donahue 2005, Briiggen et. al. 2005, Nipoti & Binney 
2005, Roychowdhury et. al. 2005, Sijacki & Springel 2006, 
and references therein). Relativistic ejecta from the AGN 
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transform into hot bubbles that eventually reach pressure 
equilibrium with the ICM and proceed to rise buoyantly 
away from the center. These bubbles are believed to be one 
of the means with which the AGN can heat the ICM. For 
example, mechanical activity associated with the inflation 
of the bubbles near the center can generate internal gravity 
waves (ie. disturbances that are subject to a restoring force 
of gravity or buoyancy) (Omma et. al. 2004) and/or sound 
waves which are believed to eventually dissipate their energy 
in the ICM. To balance cooling in a cluster of x-ray lumi- 
nosity of I/x ~ 10'''* erg s~^, a central AGN must produce 
~ 10®° erg over the entire life-time of the cluster. This is at 
the upper limit of the observed range of AGN energy output 
in galaxy clusters, based on the pV content of x-ray cavities 
(e.g. Birzan et. al. 2004). Since the pV^ energy content must 
be a lower estimate of the AGN energy output (see below), 
it is reasonable to assume that cooling is quenched by AGN 
feedback. The challenge, however, is to arrange for efficient 
energy transport from the AGN over the entire cooling core, 
or out to distances of up to ~ lOOkpc. 

Typically, studies investigating the impact of AGNs 
have tended to focus on AGN activity in only the central 
galaxy, and invoke relatively high values of thermal conduc- 
tion to redistribute the energy across the cooling core (see, 
for example, Roychowdhury et. al. 2005). Here, we propose 
instead that the cluster ICM is heated by the AGN activity 
in more than just the one central galaxy. The central re- 
gion of the clusters is rich in giant elliptical galaxies, most 
of which are presumed to harbour supermassive black holes 
that are thought to inject radio bubbles into the ICM with 
a duty cycle of ~ 10* yr per Hubble time (Best et. al. 2006). 

Upon injection, the overpressurized bubbles rapidly ex- 
pand (virtually in situ) to reach pressure equilibrium with 
their surroundings, generating waves and shocks, the dis- 
sipation of which is the primary mechanism by which the 
ICM is heated. Once in pressure equilibrium, the bubbles 
will rise buoyantly. The dissipation of the work done by vis- 
cous drag acting to retard the bubbles' motion will further 
heat the ICM, though to a much lesser extent. In response 
to the decreasing pressure, the rising bubbles will expand. 
This expansion causes a corresponding compression of the 
ICM. However, while the initial inflation of the cavity may 
trigger internal gravity waves, compression waves and shocks 
that heat the ICM, this subsequent expansion-compression 
is adiabatic in character and no heating of the ICM ensues, 
even though the ICM's internal energy will increase due to 
compression. In fact, the compression actually leads to the 
erosion of the thermal content of the ICM because the re- 
sulting increase in density results in enhanced cooling. This 
can be seen as follows. The compression of the ICM can be 
approximated as am adiabatic process so that the "entropy" 
S = T /v?^^ is conserved. The bremsstrahlung cooling time 
is tc oc T/inT"-^) ~ S°-*n-^-°^. Therefore, cooling becomes 
more efficient as n is increased during adiabatic compression. 

In our model, the heating occurs near where the bub- 
bles are produced but since the galaxies that produce these 
bubbles are distributed throughout the cluster core region, 
so too is the heating. Modest heat conduction should suffice 
to distribute further the energy over the ICM between the 



galaxies. Since the energy deposition is a local phenonemon, 
it is relatively straightforward to construct a self-regulating 
model. 

To model the interaction of the bubbles with the ICM 
(which we treat as all material not belonging to the hot and 
dilute bubbles), we use a spherical one-dimensional (ID) hy- 
drodynamical code. The code treats the bubbles and the 
ICM as a two-fluid system in a fixed dark matter halo. The 
code treats the interaction of the bubbles and ICM in a con- 
sistent semi-analytical way and assumes that the distribu- 
tion of bubbles is spherically symmetric. The code integrates 
the hydrodynamical equations of an ambient medium with 
thermodynamic properties which depend on the local filling 
factor of bubbles. 

The outline of the paper is as follows. In §2 we derive 
the equations of motion of the ambient medium and present 
details of the the interaction of the bubble fiuid with the 
ICM. In §3 we present our scheme for self-regulating feed- 
back and discuss its general properties. We describe the nu- 
merical implementation in §4 and show results in §5. We 
conclude with a summary and general discussion of alterna- 
tive heating schemes in §6. 



2 THE EQUATIONS OF MOTION FOR THE 
AMBIENT MEDIUM 

A ID hydrodynamical code cannot tackle the intrinsically 
three-dimensional (3D) problem of the evolution of individ- 
ual bubbles in the ICM, in particular if bubbles are injected 
from sources away from the center as in our proposed sce- 
nario. Therefore, a formalism for following the evolution of 
bubbles statistically must be developed. We represent bub- 
bles as a fluid specified by the following physical quantities: 
the radius of bubbles, their number density, and velocity, all 
as a function of distance from the center. The bubbles and 
the ICM are a two-fluid system that is best described as a 
single fluid which we term the ambient medium. The repre- 
sentation in terms of a single medium is applicable only if 
local pressure equilibrium between the bubbles and the ICM 
is established (Appendix A). Hereafter, quantities related to 
the ambient medium, the bubbles, and the ICM (gas outside 
bubbles) are denoted by the subscripts, a, 6, and /, respec- 
tively. Suppose that the volume filling factor of bubbles at 
r is F{r). The ambient density, p^, at position r, is defined 
as 

P. = (1 - F)p, + Fp^ , (1) 

in terms of the density inside bubbles, p^^, and the density 
of the ICM, pj , also at r. The local ambient energy, , per 
unit mass as 

y ^ (1 - F)p,u, + Fp^u^ ^ 

Pa. 

where u^^ and are the energy per unit mass of the material 
inside the bubbles and of the ICM. If the pressure, p, is 
p = if — l)wp for both the ICM and the material inside 
the bubbles. Following an analysis similar to that presented 
in Appendix A we find that the equation of state of the 
ambient medium is. 
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where 



7a -1 = 



(7i - l)(7b - 1) 



(1-F)(7,-1)+F(7, -1) 



(3) 



(4) 



and we have assumed local pressure equilibrium between 
the bubbles and the ICM. These thermodynamic relations 
supplement the following equations of motion of the ambient 
medium in a gravitational field g{r), 



dt 



1 dp 



dr 



(5) 



In the absence of dissipative heating and cooling, the adia- 
batic energy equation is, 

dw^ p dp^ 



dt 



Pa' dt 



(6) 



Dissipative heating and cooling will affect the energy equa- 
tion as described below. In the above equations the change 
of the filling factor as a result of dissipative heating and 
cooling, production of new bubbles, and bubble motion will 
have to be taken into account self-consistently. 



2.1 Motion of bubbles and their interaction with 
the ICM 

In our model, cluster galaxies affect the ICM by the pro- 
duction of over-pressurised bubbles made of hot relativistic 
plasma with an adiabatic index =4/3. Bubbles arc pro- 
duced in response to local cooling in the vicinity of galaxies 
according to the self-regulating mechanism described in §3. 

Let n^^{r), v^^{r), and R^^(r), be, respectively, the num- 
ber density of bubbles, the velocity of bubbles relative to 
the ambient medium, and the radius of bubbles, all at dis- 
tance r from the center. The motion of the bubbles is mod- 
eled as that of a pressure-free (i.e. coUisionless) fluid with 
a velocity that is determined by a balance between buoy- 
ancy and drag forces with the ICM. Bubbles are assumed to 
be injected with the same energy per bubble and the same 
initial pressure so that at distance r all bubbles have the 
same size (see equation 7) and the same velocity. There- 
fore, a hydrodynamical description for the flow of bubbles 
is self-consistent. The pressure- free assumption is also self- 
consistent since v^^ (r) is a single- valued function of r so that 
no "shocks" can form. 

The various phases of the evolution of the bubbles are 
as follows: 

• Phase I: A bubble is injected at distance r from the 
center. This over-pressurised bubble undergoes a super-sonic 
expansion in the ICM until its internal pressure drops close 
to that of the ICM. During this phase of rapid expansion, 
the bubbles heat the ICM by the production of weak shocks. 
This expansion is assumed to be quite rapid and such that 
the bubble does not move significantly away from the dis- 
tance at which it has been injected. 

• Phase 11: The dilute hot bubble then rises by buoyancy 
towards more distant regions of lower pressure. The mov- 
ing bubbles heat the ICM via viscous drag forces. As noted. 



the velocity of bubbles relative to the ambient medium is 
determined by a balance between buoyancy and drag forces 
with the ICM. Rising bubbles also expand adibatically, com- 
pressing the surrounding ICM. This process does not heat 
the ICM*. Increased ICM density does, however, elevate the 
efficiency of cooling. 

• Phase III: Finally, bubbles are destroyed by hydrody- 
namical (Raleigh- Taylor and Kelvin-Helmholtz) instabilities 
(e.g. Soker, Blanton & Sarazin 2002). 

2.1.1 Initial parameters of an injected bubble 

Let a bubble be injected with an initial energy E^^- and an 
initial internal pressure p^^.. Since E^^. = u^^-p^^ AnR^^/^/3, 
where Pj, . and - are the initial mass density and energy 
per unit mass in the bubble, the equation of state p^, = 
(7b - l)PbWb yields. 



-(7 



1)^ 

Pbi 



1/3 



(7) 



2.1.2 Bubble radius as a function of distance 

The bubble radius, Ri,{r), at any distance r from the center 
can be estimated from the adiabatic condition PbPb"^'' = 
const. This condition must hold in Phase I and Phase II 
because the speed of sound of the relativistic plasma in the 
bubble is much larger than its expansion velocity. The adi- 
abatic condition yields 



^« = -Rbi 



Pb 



p{r) 



3Tb 



(8) 



which for — 4/3 gives the weak dependence, R^^ oc 
{Phi/Py^*' on the pressure ratio. 



2.1.3 Heating by shocks 

Wo estimate the energy transferred from bubbles to the 
ICM by means of weak shocks during Phase I. The equa- 
tion of state of the material inside the bubble gives E^^ = 
(47r/3)-R^,'Vb/(7b ^1) for the thermal energy of a bubble of 
radius R-^ with internal pressure p,^ . At the end of the rapid 
expansion in Phase I we assume that p^^ ~ p(r) where p(r) 
is the ICM pressure at the distance at which the bubble has 
been injected. Substituting p = p(r) in (8) to get R^^ , we 
find that the thermal energy at the end of Phase I is 

This equation implies that E^^ oc pRj^ can be substantially 
smaller than the initial energy of the bubble. The heat trans- 
ferred to the ICM by the generation of shocks is 



* Deviation from adiabatic compression may result in bulk mo- 
tions in the ICM. These motions can eventually dissipate into 
heat via turbulence or ordinary viscosity. The maximum heating 
rate that is obtained this way is slightly less that the drag heating 
rate discussed in §2.1.4 
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.(10) 



where < C^^^ < 1 is an efficiency parameter. In the present 
treatment, wo assume that the energy transferred to the 
ICM is deposited iocaliy. The expression (10) does not take 
into account the energy needed to place the bubble with 
initial radius R^^. and energy E^^. into the ambient medium. 
For a given E^^. , this energy is negligible if i?^; is sufficiently 
small. 



2.1.4 Buoyancy and heating of the ICM by drag with 
rising bubbles 

The dynamical evolution of the system comprised of an 
expanding, rising bubble in a bubbly medium subject to 
radiative cooling is, in general, a highly non-trivial prob- 
lem, involving a host of complex phenomena including the 
compressible nature of the bubbly medium, the evolution 
and dissipation of shocks, wakes and waves in this medium, 
bubble-bubble interactions mediated by bubble wakes, bub- 
ble coalescence, changing bubble shapes and orientations, 
etc. We write the equations of motion of a bubble in Phase 
II. We restrict ourselves to a simple situation of a system 
involving a single spherical bubble rising along the radial 
direction in the ICM. The equation of motion governing the 
bubble's rise is 
d_ 
dt 



(M„«J=F,„ 



(11) 



where the two terms on the right represent the buoyancy 
force acting on the bubble of radius and the drag force on 
the bubble as it moves through the ICM. The left-hand term 
is the rate of change of the Kelvin impulse (also known as 
the added-mass force) corresponding to the displacement of 
the ambient medium due to the bubble's rising motion. Ana- 
lytic expressions for the "virtual added-mass" of the bubble 
can be obtained for simplified motions. For a laminar flow 
generated by a spherical bubble moving in an incompress- 
ible fluid the expression for the virtual added mass is (e.g. 
Kendoush 2003), 



Mo 



Cam—R^^P^, Cam 



(12) 



where R is the rate of expansion of the bubble. We have 
neglected the inertia of the bubble as it is much smaller 
than the added-mass term. 

A bubble of radius R^ at distance r experiences a buoy- 
ant force of 



47r 3 



(13) 



where we have assumed that ^ Pa- Here, is the grav- 
itational force field (per unit mass) as measured in a frame 
of reference attached to the ambient medium at distance r. 
The actual gravitational field, g, differs from g^^^ only when 
the system deviates from hydrostatic equilibrium, such as in 
a cooling runaway where g^^^ « 0. We write the drag term 
on the bubble as it moves through the ICM with velocity v-^^ 
as 



F,.., = ^C^R^'P.V^^ ■ (14) 

where C^ is the drag coefficient and we have assumed that 
<SC Cs, In our model all bubbles at the same distance share 
the same velocity. But in a realistic situation, bubbles are 
generated with different initial energies and at random direc- 
tions. So a bubble may collide with others. Hence, to model 
this effect we write instead of P; in the drag equation. 

Typically, the magnitude of added-mass term is much 
smaller than the buoyancy and the drag terms and there- 
fore, to first order, the rise velocity of the bubbles can be 
estimated by equating the drag and buoyancy forces (e.g. 
Churazov et. al. 2001), yielding: 



2 _ ^g^aRb 



3a 



(15) 



The bubble radius, Ri,{r), as given by (8) depends only on 
p{r). Therefore, all bubbles with the same initial radius and 
internal pressure share the same velocity at r, independently 
of the distance at which they have been injected into the 
ICM. When is non-negligible compared to the speed of 
sound in the ambient medium, the drag becomes more effi- 
cient. Here we account for the increase in the drag coefficient 
by using the following expression for 



3a 



1 + 



8fleff-^b 

sac? 



(16) 



As a bubble rises, the gravitational potential energy of 
the system (the bubble and the surrounding medium) de- 
creases. A Rising bubble also expands and its internal en- 
ergy decreases. The total energy available as the bubble rises 
from position 1 (point of injection) to position 2 is 

Y^b'pa5eff dr . (17) 



where the first difference pair, corresponding to the change 
in the bubble's internal energy, it can be evaluated using 
equation the condition p/p^l = const, for adiabatic changes 
of the bubble fluid as 



{Ehi — E\j2) = -Ebi ( — 



1-1/7. 



1- - 



l-l/7h 



(18) 



Assuming that the ambient medium is in hydrostatic equi- 
librium, the second term too can be straightfowardly evalu- 
ated, yielding: 



(19) 



where is the gravitational potential. Consequently, AE2 = 
(1 + 7b)(^bi ~ Eh2). This is the energy available to lift the 
bubble against the drag force, overcome the inertia of dis- 
placing the surrounding fluid, as well as compress the ambi- 
ent medium. 

In general, it is not possible to determine how much of 
this energy that will go towards compressing the ambient 
medium and how much will towards heating the ICM, be it 
due to drag or through the eventual dissipation of the walces 
and motions associated with the displacement of the ICM, 
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unless one knows the details of the rise-expansion process. In 
the limit that the bubble rises slowly, maintaining pressure 
equilibrium with its surroundings at all times, the bubble ex- 
pansion/ICM compression is adiabatic and the change in the 
internal energy of the bubble goes entirely towards raising 
the internal energy of the ambient medium. In this case, the 
total energy available for heating the ICM is simply that 
due to the change in the potential energy and the corre- 
sponding heating rate (in units of energy per unit time) is 
'^b-^bu ~ (47i"/3)-Rb ^Po-^off '^b • To get the heating rate per 
unit mass, 'H^^^^^ , we multiply this expression by the num- 
ber density of bubbles, n^^ , and divide it by the mass density 
of the ambient medium per unit volume (see equation 1). 
The result is 
47r . 



(20) 



It is instructive to relate this heating term to the physi- 
cal parameters of the bubbles and their injection rate. As- 
sume that bubbles are injected within a sphere of radius R.^. 
around the center, at a rate N-^. per unit volume. Neglecting 
bubble destruction, the total number of bubbles within a dis- 
tance r » Rj^j is -^Vj„j(47r/3)-Rj„j^r/«jj, so that the number 

density within r is 7i^{< r) ^ -^inj^inj^/(^b''^)- Approxi- 
mating n^(r) ~ n^{r) in (20) and using (15) gives 



(21) 



It is interesting that this does not depend on the drag coeffi- 
cient, Cj. We can further simplify this expression by noting 
the bubble energy at r is -Bj^ AnRf^^p/3/{'y^ - 1). The final 
result is 

Ko. -(7b-l)iV„,5.«%^^. (22) 

For r < R.^. the heating rate is obtained by replacing R^^. 
with r in the last expression. 

In the rmmcrical code, the heating is computed explic- 
itly using (20) . We do so because apart from the initial accel- 
eration phase, the bubble's motion should be well described 
by the "slow rise" approximation. In the worst case scenario, 
we are underestimating heating rate in Phase Ilhy a factor 
of 1 + l/7j^. This uncertainty is not of concern because in 
our adopted scenario where the bubbles are initially injected 
in a highly overpressurized state (pi <C Pbi); the heating of 
the ICM during Phase /dominates over the heating by the 
rising, expanding bubble. 



3 SELF-REGULATING FEEDBACK 

We write the energy equation in terms of the ICM entropy 



5; = Uj/pj'"'! as 



Pi 



Tl- 



'5, = ?i 



^cool 



(23) 



where icooi is given in Appendix B. As explained before, 
two sources contribute to H: a) heating by weak shocks 
generated during the initial rapid expansion of the over- 
pressurized bubbles, and b) heating by drag with the ICM 



as bubbles rise up by buoyancy. Adiabatic work done by 
the bubbles as they move upward does not change <Sj and 
therefore does not contribute to H. We introduce R.^. , the 
distance from the center within which bubbles arc injected. 
At r > R-^. , bubbles can be present only as a result of flow 
of bubbles from regions with r < R^^. . Given the injection 
rate (number per unit volume per unit time), iV(^j(f, r), of 
bubbles, we write the local heating rate (energy per unit 
mass of the ICM) via weak shocks as 



H = N 

' "-wsh ' inj 



(24) 



where AEi is given by (10). As a recipe for self-regulated 
feedback, we assume that bubbles are generated only if 
is negative. We write the flux of injected bubbles at a point 
r < R,„. as 



N. . = 



V 'Si 2 



for 5, < 



(25) 



and zero, otherwise. The free parameter rj represents the 
product of the mass fraction of cold gas that accretes onto 

the AGNs and the efficiency of AGNs at transforming the 
accreted mass into bubbies. The cofd gas is assumed to be 
generated at a rate (in units of mass per unit volume per 
unit time) of p^S^/S^. This cold gas is assumed to form by 
condensation via the development of thermal instabilities 
(e.g. Field 1965). Finally, c is the speed of light and, as 
before, iSj^j is the energy with which a bubble is injected. 

Wc can use the above formalism to investigate the role 
of heating by weak shocks, the dominant source of heating 
in our model. Substitute ASi and N-^^-^ from (10) and (25), 
respectively, into (24) yields 



TL = —Uj^ ^ , where tt^ = C^^^r](? 



1-1^ 



,(26) 



where p = p{r) is the ambient pressure at r and is the 
internal bubble pressure with which the bubbles are injected. 
This expression holds for iSi < 0. For Sj > 0, we have H = 0. 
Substituting (26) in (23) and solving for Sj we get. 



'?~cool 



where 



'^cool 



fl+^ 

V u 



(27) 



(28) 



The above equation implies that cooling is substantially sup- 
pressed at low temperatures, w <S i while at higher tem- 
peratures, shocks are unable to suppress cooling. For an ICM 
at temperature corresponding to velocity V an efficient sup- 
pression of cooling at distance r occurs if 



(29) 



where we have assumed ^ p. Note that since N^^. oc 
neither the shock nor the drag heating depend on 
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4 NUMERICAL MODEL 

We develop a numerical hydrodynamical model to describe 
the evolution of the ICM and the bubbles in a fixed dark 
matter halo. The density profile in the halo is approximated 
by the NFW profile (e.g. Hayashi et. al. 2004) of the form 

where s = r/Rv and Rv is the virial radius of the halo. 
The parameter 5c is determined by the condition that the 
mean density within r = i?v is 178 times the background 
density, p{t) . We adopt the value c = 5 for the concentration 
parameter (e.g. Balogh et. al. 2005). 

The ambient medium is represented by lagrangian shells 
of fixed mass. In each time-step, the equations of motion of 
the bubble and the ambient fiuids are solved in the following 
order: 

(i) advance the ambient shells according to equations (5) 

and (6) without radiative cooling or dissipative heating. 

(ii) update bubble properties in each shell as a result of the 
change in incurred in the ambient medium (see Appendix 
A) after the hydrodynamic step (i). 

(iii) estimate dissipative heating (by weak shocks and drag) 

and radiative cooling of ICM and update the corresponding 
ambient pressure and bubble properties without moving the 
ambient shells. 

(iv) inject the ambient medium with new bubbles according 
to the recipe outlined in (§3). 

The integration of the hydrodynamic equations (5) and 
(6) includes artificial viscosity (e.g. Richtmycr & Morton 
1967; Thoul & Weinberg 1995; Nusser & Pointecouteau 
2006) in order to deal with any abrupt changes in the veloc- 
ity which may arise if a cooling runaway develops. The bub- 
ble fluid is moved through the shells of the ambient medium 
using a first order flux conserving scheme. The velocity of 
this fluid relative to the shells is given by (16). Finally, a 
leapfrog time integration scheme is used. 



mass ratio at r = Rv is taken to be the background ratio 
nb/fim. The mass in each ambient shell is then estimated 
from the gas density profile. Because of the limited num- 
ber of shells, a constant temperature does not necessarily 
guarantee that the shells are actually in strict hydrostatic 
equilibrium. Therefore, given the mass in each shell, we use 
the discrete form of the equations of motion as they ap- 
pear in the numerical hydrodynamical code (e.g. Weinberg 
& Thoul 1995) to obtain a temperature (thermal energy per 
baryon) in each shell such that strict equilibrium is achieved 
at the initial time. This yields nearly isothermal tempera- 
ture profiles with deviations that increase with radius. 

In all runs, the code evolves the shells from redshift 
z — 2 until z = 0. In some runs with insufficient feedback to 
counter cooling, the hydrodynamic time-step becomes ex- 
ceedingly small. Whenever this happens, the shells colder 
than 1 keV and denser than 1.5cm~^ are held fixed and are 
not moved by the code at any later time. 

Results from simulations for clusters with virial radii 
of _Rv = 4 Mpc and 7?v = 2 Mpc are shown in Fig. 1 and 
Fig. 2, respectively. In the column to the left we show results 
which include radiative cooling and no heating. For both 
values of _Rv, effects of cooling become significant already 
byz = l. Atz<0.5a cooling runaway develops in the 
inner regions: cooling becomes increasingly more rapid due 
to the enhanced density and lower temperature. 

The effect of heating by bubbles, according to our feed- 
back recipe described in §3, is explored in the middle and 
right panels for two values of the efficiency parameter r). To 
model the effect of multiple AGN activity, bubbles are al- 
lowed to be injected anywhere in the ICM (i.e. R^^^ = R^). 
As noted before the heating rate of the ICM does not de- 
pend on the initial bubble energy, E^^^ . All results in these 
two columns are for C . = 1 and C =1, and an initial bub- 

d wsh ' 

ble pressure of 10^ times the central ICM pressure. Cooling 

fiows for the two different clusters arc completely suppressed 
for 77 as low as 0.01, while a value r] — 10^'' gives moderate 
suppression of cooling as is required by the observations. 



5 RESULTS 

We adopt a flat cosmology with baryonic and dark matter 
density parameters respectively Ob = 0.044 and Om = 0.27, 
a cosmological constant corresponding to a density param- 
eter of Ha = 0.686, and a Hubble constant of Ho — 
70kms-^Mpc-^ (e.g. Spergel et. al. 2003). 

The number of shells describing the ambient medium is 
50 in all of the runs. The initial shells are placed at loga- 
rithmically spaced distances from the center with the near- 
est shell r — 15kpc and the farthest lying at the virial ra- 
dius of the halo. The initial temperature proflle is assumed 
to be isothermal with a temperature, T (in keV), given by 
1.16 X 10'^ ksT / finip = (2/3)GMv/7?v where fce is the Boltz- 
marm constant, fi is the mean molecular weight and nip is 
the mass of the proton. The equation for hydrostatic equilib- 
rium implies a slope of —9/4 for the gas density at r = R^. 
The gas density at any r {r < Rv) is solved numerically us- 
ing the equation of hydrostatic equilibrium including gas 
gravity. In solving this equation, the gas to dark matter 



6 SUMMARY AND DISCUSSION 

We propose that modest radio outbursts from many cluster 
galaxies may help to quench cluster cooling flows. We ap- 
peal to typical massive elliptical galaxies which arc thought 
to have radio outbursts with a duty cycle of ^ 10** yr per 
Hubble time (Best et. al. 2006) and are concentrated in the 
cluster core (within the central 200 kpc or of order twice 
the extent of the cooling region). Weak shocks generated 
during the initial bubble expansion and viscous/turbulent 
drag on buoyantly rising bubbles are the principal dissipa- 
tive heating processes. Indirect evidence for multiple heating 
sources in the form of AGN comes from the frequent pres- 
ence of x-ray "ghost" cavities that survive long after any 
associated radio lobes have decayed, and from amorphous 
radio halos that further suggest an ICM percolated with a 
distribution of bubbles whose contrast in the X-ray heis been 
diminished due to adiabatic expansion (cf. Heinz & Chura- 
zov 2005). Our scenario is in agreement with recent obser- 
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Figure 1. The ICM electron number density (top), temperature (middle) and entropy (bottom) profiles versus distance from the cluster 
center (in Mpc). The left column shows profiles obtained with radiative cooling and no energetic feedback. In the middle and right panels 
feedback is done via bubbles generated in the entire cluster with efficiency parameter r/ = lO^'^ and r; = 10^''', respectively. The solid, 
dotted, dashed and dash-dotted curves correspond to redsliifts z = 1.6, 1, 0.5, and 0, respectively. The results correspond to an NFW 
dark halo profile with a virial radius of = 4 Mpc and a concentration parameter of c = 5. The number density, temperature, and 
entropy are given in cm~^, keV, and keVcm^, respectively. 
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vational evidence for the existence of AGN activity in about 
5% of cluster galaxies at low redshift (Martini et. al. 2006). 
Cosmological evolution of the AGN inevitably increases the 
feedback at earlier epochs (2: ^ 2). Our scenario alleviates 
some of the problems that may arise if energy feedback is 
generated by a single central source (AGN) harboured in 
the central cD galaxy. The needed energy transport over 
the entire cooling flow region is not easy to provide if the 
energy source is associated only with the central AGN (e.g. 
Mathews, Faltcnbachor & Brighcnti 2005). One approach is 
to ensure that mechanical disturbances are able to prop- 
agate large distances before they are damped by viscosity 
(e.g. Ruszkowski, Briiggen & Begelman 2004; Reynolds et. 
al. 2005). This approach has two potential shortcomings. 
Firstly, the heating rate per baryon of dissipating waves 
should decay with distance from the central source as . 
This is too steep a gradient for the heating to affect the ICM 
over the entire region that is susceptible to cooling. In fact, 
Roychowdhury et. al. (2004) find that the heating leads to a 
convectively unstable core profile that would potentially de- 
stroy any pre-existing metallicity gradient, in conflict with 
the x-ray observations, and to remedy this, Roychowdhury 
et al. (2005) invoke relatively high values of thermal con- 
duction to effect rapid transport of the energy out of the 
center. Secondly, the viscosity has to be finely tuned: too 
high a value results in dissipation much to closely concen- 
trated near the central source, and too low a value means 
the heating is inefhciont. Since the relevant viscosity coeffi- 
cient scales as T^^'^ (Spitzer 1962), it is difficult to envison 
that a viscous solution that relies on fine-tuning prevails in 
all clusters. Moreover, the ICM is weakly magnetized, and 
the nature of transport processes, like thermal conduction 
and viscosity, in such media is not well understood. 

The model presented here heats the ICM through 
shocks generated during the initial inflation of the bubbles, 
and drag forces acting on the rising bubbles. However, the 
latter relies on some viscosity to dissipate heat but it is 
not an important source of heating and with respect to the 
former, we do not need to fine-tune the viscosity since the 
energy distribution is effected by multiple AGNs hosted in 
a spatially extended distribution of galaxies. Moreover, it 
may be possible to avoid the issue altogether. Heinz & Chu- 
razov (2005) have recently argued that in a bubble-filled 
medium of the kind proposed here, the bubbles themselves 
act as catalysts to convert shock and waves into heat. Addi- 
tionally, the interaction between the bubbles and the strong 
shocks will both broaden and distort the shock fronts, caus- 
ing the features to appear weaker in X-ray observations and 
resulting in mistaken assumptions that cavity inflation is a 
gentle process and that the ICM is not subject to strong 
shocks. Both of these features further contribute to making 
our model highly viable. 

An additional principal feature of our model is the in- 
corporation of self-regulation between heating and cooling. 
This should inevitably lead to a correlation between the me- 
chanical heating luminosity and the cooling rate. The former 
is measured by the properties of x-ray cavities and the latter 
by the cluster x-ray luminosity. Evidence for such a correla- 
tion is given by Birzan et. al. (2004). 



In fact, the estimated ratio of mechanical to x-ray lumi- 
nosity in this latter paper falls short of that needed by Cro- 
ton et. al. (2005) in order to quench cooling flows on group 
and cluster scales by a factor of 3-10 (Best et. al. 2006). 
However, the inferred mechanical luminosities are based on 
cavity pV estimates. In our model, a primary source of en- 
ergy transfer to the ICM is via weak shocks generated by 
the expansion of the initially over-pressurised bubble. The 
initial bubble energy is related to the cavity pV by equa- 
tion (9), and the resulting injected mechanical energy may 
exceed pV by a factor of up to 10. Note that this equation 
neglects the work done by the AGN jet to place the initial 
seed bubble in the ICM, which can be dissipative and con- 
tribute to heating the ICM if the seed bubble is ejected with 
supersonic speeds. Furthermore, these estimates neglects the 
energy content in waves detected in a few clusters (Fabian 
et. al. 2003; Nulsen et. al. 2005, McNamara et. al. 2005). 
Our model consequently helps in resolving the apparent dis- 
crepancy in the estimated feedback required to halt cooling 
flows even in massive clusters. 

An alternative heating source involves cosmic rays in 
the ICM cooling core, injected by a central AGN (Rephaeli 
and Silk 1995). The present model provides a distributed 
source of cosmic rays whose observable signatures could also 
include a Sunyaev-Zeldovich signal (Pfrommer & Ennslin 
2004) as well as an additional hadronic heating source (Co- 
lafrancesco, Dar & De Rujula 2004) that might lead to a de- 
tectable gamma ray flux. More exotic heating sources invoke 
collisions between superheavy dark matter particles and pro- 
tons (Qin & Wu 2001, Chuzhoy & Nusser 2004) as well as 
by the products of dark matter neutralino annihilations (Co- 
lafrancesco 2004). 
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APPENDIX A: THERMODYNAMICS OF THE 

BUBBLES AND ICM 

Al Numerical description of bubbles in ambient 
medium 

For simplicity the mass of bubbles is neglected. 

As bubbles are transported in the ICM in an infinites- 
imal time step, the ICM in a shell is either compressed or 
diluted depending on whether or not number of bubbles in 



the shell has increased or decreased. This is a purely adia- 
batic process that conserved the entropies inside a bubble 
and in the ICM. 

Consider a shell with pressure pi containing bubbles 
each of radius Rhi. Assuming that 5N additional bububles 
enter the shell, we would like to estimate the ambient pres- 
sure, p, the radius, Rh, and the filling factor, F, after the 
bubbles enter the shell. For concreteness let 5N > 0, i.e., 
there is a net increase in the number of bubbles in the shell. 
Before these bubbles enter the shell, the total material in 
the shell occupies a volume Vi = 47r(r?|_i — rf)/3 of which 
Vbi = FiVi is taken by the bubbles and Vn = {1 - F)Vi by 
the ICM. 

After 5N bubbles enter the shell, the same original ma- 
terial occupies a volume Vi - 5V where 5V = 5N4nRj/3 
with JZ^ being the new radius of a bubble in the shell. 
The entropies in the bubbles and the ICM of the original 
arc conserved so that pVj^ = PiVjl and pV^^ = PiVti^t 
where Vi and Vb are the volumes occupied by the ICM and 
bubbles of the original material after the 5N bubbles en- 
tered the shell. Using these entropy conservation relations 
in Vi — 6V = Vi+Vh gives 



211 '+(!■ 

p 



(Al) 



Entropy conservation inside a bubble gives 7?^^ = 
i?t, j(pi/p)^/Tb which after substituting in 6V = dN4nR^,^/3 
yields. 



l=fF,+.iV±^7Z,?)f^)^%(l-F,)(^) ■ .(A2) 



The bubble radius and their filling factor after the change 
can easily be computed from the adiabatic conditions once 
the pressure, p, is known. 



Al.l Effect of heating and cooling on the ICM and bubbles 

As mentioned previously, the ICM is meant to describe the 
inter-cluster plasma lying outside the bubbles. As an in- 
finitesimal amount of energy is removed from the ICM by 
cooling or added, we must describe the reaction of the ICM 
and bubbles to such infinitesimal change of the energy. Let 
an energy 5E be added as heat to the ICM in the shell 
i + 1/2. For concreteness we assume SE > 0, but the same 
reasoning applies to cooling SE < 0. In the calculation we 
assume that the that the system reacts in the following two 
step sequence. First the pressure of the ICM is changed to 
Pi + (7i - ^)SE/Vo where Vo = 47r(rf+i - rf )/3 before the 
bubbles are compressed. Next, the under- pressured bubbles 
are compressed by the ICM. The expansion of the ICM and 
the compression of the bubbles in the second step are both 
assumed to be adiabatic processes. The following adiabatic 
conditions hold in the second step. 



Pi + (7i - 1) 



SE 



and 



pVb^ =Pin> (A3) 
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where p, Vi, and 14 are the pressure, the volume of the 
ICM, and the volume of bubbles at the end of the second 
step. These conditions give the following equation for p 

' .(A4) 



APPENDIX B: THE COOLING FUNCTION 
AND THE X-RAY EMISSIVITY 

Wo approximate the cooling time by (Balogh, Babul & Pat- 
ton 1999) 

'"°°'-p,[l + B2^/T,(r)] ' ^^'^ 

where, ji is the mean molecular weight, Bi — 3.88 x 
10"sK-i/2cm-3, B2 = 5 X and /m is a metallicity- 

dependent constant that is 1 for solar metallicity and 0.03 
for zero metallicity. 

The volume X-ray emissivity of the ambient medium is 

2 /xmn tcooi 
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